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ABSTRACT 



Context. The red hypergiant NML Cyg has been assumed to be at part of the Cyg 0B2 association, although its distance has never 
been measured directly A reliable distance is crucial to study the properties of this prominent star. For example, its luminosity, and 
hence its position on the H-R diagram, is critical information to determine its evolutionary status. In addition, a detection of the radio 
photosphere would be complementary to other methods of determining the stellar size. 

Aims. We aim to understand the characteristics of NML Cyg with direct measurements of its absolute position, distance, kinematics, 
and size. 

Methods. We observe circumstellar 22 GHz H2O and 43 GHz SiO masers with the Very Long Baseline Array to determine the 
parallax and proper motion of NML Cyg. We observe continuum emission at 43 GHz from the radio photosphere of NML Cyg with 
the Very Large Array. 

Results. We measure the annual parallax of NML Cyg to be 0.620 + 0.047 mas, corresponding to a distance of 1.61^^ }] kpc. The 
measured proper motion of NML Cyg is yU;^ = -1.55 + 0.42 mas yr~^ eastward and fiy = -4.59 + 0.41 mas yr~^ northward. Both the 
distance and proper motion are consistent with that of Cyg 0B2, within their joint uncertainty, confirming their association. Taking 
into consideration molecular absorption signatures seen toward NML Cyg, we suggest that NML Cyg lies on the far side of the 
Cyg 0B2 association. The stellar luminosity revised with our distance brings NML Cyg significantly below the empirical luminosity 
limit for a red supergiant. We partially resolve the radio photosphere of NML Cyg at 43 GHz and find its diameter is about 44 
mas, suggesting an optical stellar diameter of 22 mas, if the size of radio photosphere is 2 times the optical photosphere. Based on 
the position of circumstellar SiO masers relative to the radio photosphere, we estimate the absolute position of NML Cyg at epoch 
2008.868 to be Qfj2ooo = 20M6"^25?5382 + 0?0010, 6nooo = 40°06'59:'379 + 0:'015. The pecuhar motions of NML Cyg, the average 
of stars in Cyg 0B2, and four other star-forming regions rules out that an expanding "Stromgren sphere" centered on Cyg 0B2 is 
responsible for the kinematics of the Cygnus X region. 

Key words, astrometry — masers — parallaxes — proper motions — stars: individual (NML Cyg) — supergiants 



1. Introduction 



NML Cyg is one of the most massive and luminous red 
hypergiants. It is a sem i-regular sta r with period ^ 1000 
days ( Denisenkova et al.l J 988; Monnier et all Il997h and has 
been suggested to be associated with Cy g QB2, possi b ly the 
largest stellar association in the Galaxy ("Habing e t all 1 19821: 
Im orris & Jural |l983). Previous distance estimates for NML 
Cyg have been mainly based on the assumption that it is at 
the distance with Cyg OB 2. Even so, the distance to Cyg 
OB 2 is hard to determine accurately by photometry and spec- 
troscopy, due to difficulty in calibrating spectral types and lu- 
minosities. Previous investigators have derived distance moduli 
of 10.8 - 11.6 (IJohnson & Morgan! 11954: Redd ish etal.l 11966 : 
[Walbornl "197 3"; 'Abbott et al.','l98l';"Torres-Dod2en et al.l 11991 
Massey & Thompson, 1991; Hanson, 2003; Kharchenko et al 
1200 5), corresponding to distances of 1.45 - 2.10 kpc to Cyg 
0B2. These values are margi nally differen t from 1.22 ± 0.3 kpc 
to NML Cyg determined bv IPanchi et al.l (l2Q01i) . based on the 
comparison of Doppler velocities and proper motions of two dis- 



crete dust shells around NML Cyg measured with interferometry 
at 1 1 yum over a period about 6 yr. 



Adopting a d is tance of 1.74 + 0.2 kpc 
(iMassev & ThompsonL Il99lh . ISchuster et aP (l2QQ9l) esti- 
mated NML Cyg's minimum bolometric luminosity to be 
3.15 ± 0.74 X 10^ Lo, which is similar to that of other 
red hypergiants and places it near the empirical upper- 
l uminosity boundary i n the Hertzsprung-Russell (H-R) diagram 
(ISchuster et al.L l2006l) . However, the estimated luminosity of 
a star depends on the square of its distance. Thus, an accurate 
distance is crucial to derive a reliable luminosity and the 
location on the H-R diagram. Recently, trigonometric parallaxes 
of circumstellar maser sources (H2O or SiO) surrounding red 
hypergiants with VLBI phase-r eferencing have been measured 
with accuraci es of 20 - 80 yuas (IChoi et al.L l2008l: IZha ng et aP. 
2012; Asaki e t al.L 1201 0). We carried out a program to measure 
the trigonometric parallax and proper motion of masers in the 
circumstellar envelope (CSE) of NML Cyg with the National 
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Radio Astronomical Observatory'^ (NRAO's) Very Long 
Baseline Array (VLB A). 

In addition to distance, another fundamental stellar param- 
eter is size. NML Cyg's high mass-loss rate results in a dense 
circumstellar envelope, and the star is faint and hard to observe 
at visual wavelengths due to high e xtinction. However , it is ex- 
tremely luminous in the infrared. Ridgway et al.' ('1986) used in- 
frared speckle interferometric observations of NML Cyg and a 
simple model of multiple shells to derive an inner radius for 
its dust shell of 45 ± 12 mas. They also obtained the star's ef- 
fective temp erature, Jeff ^ 3250 K, which is app ropriate for an 
M6 III star (John sonL[l967h . lMonnier et al.l(ll997l ) conducted in- 
frared interferometry at 11.15 /im, yielding evidence for mul- 
tiple dust shells and asymmetric dust emission around NML 
Cyg. They obtained better fits to all data by assumin g a Te ff 
^ 2500 K, u s ing a model (iRowan-Robinson & Harrisl Il983\ 
[Blocker et al.l (l200lh presented diff'raction-limited 2.13 fim ob- 
servations with 73 mas resolution and obtained a bolometric flux 
of Fhoi = 3.63 X 10^ W m"^, corresponding to a stellar luminosity 
of Lboi = 1.13 X 10^- (J/kpc)^ Lq. Using the spectral energy dis- 
tribution (SED) from 2 to 50 yum, they derived a stellar diameter 
of 16.2 mas assuming Tq^ = 2500 K. Clearly, a direct detection 
of the star would be complementary to the methods mentioned 
above. Radio continuum emission from the evolved star's pho- 
tosphere can be imaged with the Very Large Arra y (VL A) , using 
H2O or SiO masers as a phase reference (Reid & Menten [ ll990l 
[1997I I2OO7I: IZhang et al.[ 12012 ). Such observations aUow us to 
directly measure the size of NML Cyg's radio photosphere and, 
using the observed flux density and our measured distance, its 
efl'ective temperature. 
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Fig. 1. Images of the H2O maser reference spot and the extra- 
galactic radio sources used for the parallax measurements of 
NML Cyg at the middle epoch (2008 November 08). Source 
names are in the upper left corner and the restoring beam (dot- 
ted ellipse) is in the lower left corner of each panel. Contour 
levels are spaced linearly at 3 Jy beam"^ for NML Cyg and 1.5 
mJy beam"^ for the background sources. 



2. Observations and data reduction 

2.1. VLB A phase-referencing observations 

We conducted VLBI phase-referencing observations of NML 
Cyg at 22 and 43 GHz with the VLBA under program BZ036 
with 7-hour tracks on 2008 May 25, August 13 and November 
8, and 2009 January 27 and May 3. This sequence samples the 
peaks of the sinusoidal trigonometric parallax curve and has low 
correlation coefficients among the parallax and proper motion 
parameters. We observed several extragalactic radio sources as 
potential background references for parallax solutions. We alter- 
nated between two 16 min blocks at 22 and 43 GHz. Within 
a block, the observing sequence was NML Cyg, J2052-h3635, 
NML Cyg, J2044-h4005, NML Cyg, J2046-h4106, NML Cyg, 
J2049+4118. We switched between the maser target and back- 
ground sources every 40 s for 22 GHz and 30 s for 43 GHz, 
typically achieving 30 s for 22 GHz and 20 s for 43 GHz of on- 
source data. 

We placed observations of two strong sources 3C345 
(J1642-h3948) and 3C454.3 (J2253-hl608) near the beginning, 
middle, and end of the observations in order to monitor delay and 
electronic phase differences among the intermediate-frequency 
bands. The rapid- switching observations employed four adja- 
cent bands of 8 MHz bandwidth and recorded both right and 
left circularly polarized signals. The H2O and SiO masers were 
contained in the second band centered at a Local Standard of 
Rest (LSR) velocity of -7 km s"^ In order to do atmospheric 
delay calibration, we placed "geodetic" blocks before and af- 
ter our phase-reference observations. These data were taken in 

^ The national Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc 



left circular polarization with eight 8 MHz bands that spanned 
480 MHz of bandwidth between 22.0 and 22.5 GHz; the bands 
were spaced in a "minimum redundancy configuration" to uni- 
formly sample, as best as possible, all frequency diff'erences. 

The data were correlated in two passes with the VLBA cor- 
relator in Socorro, NM. One pass generated 16 spectral channels 
for all the data and a second pass generated 256 spectral chan- 
nels, but only for the single (dual-polarized) frequency band con- 
taining the maser signals, giving a velocity resolution of 0.42 
km s"^ for 22 GHz and 0.21 km s"^ for 43 GHz, respectively. 
Generally, the data calibration was performed with the NRAO 
Astronomical Image Processing System (AIPS), f ollowing sim- 
ilar procedures as described in Reid et al] (12009a'). We used an 
H2O maser spot at a Vlsr of 5.22 km s"^ and a SiO maser spot at 
a Vlsr of -5.70 km s"^ as phase-reference sources, because they 
are considerable stronger than the background source and could 
be detected on individual baselines in the available on- source 
time. We imaged the calibrated data with the AIPS task IMAGR; 
Fig. [T] shows example images at 22 GHz from the middle epoch. 
Table [T] lists the positions, intensities, source separations, LSR 
velocity of the reference maser spot and synthesized beam sizes. 

2.2. VLA observation for radio photosphere 

In addition to the VLBA observations, we carried out an observa- 
tion of the radio photosphere and maser emission toward NML 
Cyg at 43 GHz with the VLA in its largest (A) configuration 
under program AZ178 on 2008 December 20. We used a dual 
intermediate frequency band setup in continuum mode with a 
narrow (3. 125 MHz) band centered at Vlsr of - 1 8 km s"^ for the 
V = 1, / = 1 ^ SiO maser line (rest frequency of 43122.027 
MHz) and a broad (50 MHz) band centered ^50 MHz above in 
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Table 1. Positions and Brightnesses 



Source 


R.A. (J2000) 


Dec. (J2000) 




^sep 


P.A. 




Beam 






I ) 


\jy/ ucdiii J 


\ ) 


\ ) 


(km s~^) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


22 GHz 
















NML Cyg 


20 46 25.5444 


40 06 59.383 


16-77 






5.22 


0.6 X 0.3 @ -08 


J2052+3635.... 


20 52 52.0550 


36 35 35.300 


0.050 


3.7 


-19 




2.7 X 1.4 @ +31 


J2044+4005.... 


20 44 n.0877 


40 05 36.360 


0.005 


0.4 


+87 




2.7 X 1.3 @ +25 


J2046+4106.... 


20 46 2L8414 


41 06 01.107 


0.014 


1.0 


-01 




2.7 X 1.2 @ +24 


J2049+4118.... 


20 49 43.9548 


41 18 14.583 


0.011 


1.3 


+28 




2.7 X 1.2 @ +26 


43 GHz 
















NML Cyg 


20 46 25.5378 


40 06 59.413 


3.4 






-5.70 


0.3 X 0.1 @ -07 


J2046+4106.... 


20 46 2L8414 


41 06 01.107 


0.017 


1.0 


-01 




0.7 X 0.4 @ -25 



Notes. The fourth and seventh columns give the peak brightnesses {S p) and Vlsr of reference maser spot. The fifth and sixth columns give the 
separations {Osep) and position angles (P.A.) east of north of the background sources relative to the maser. The last column give s the FWHM size 
and P.A. of the Gaussian restoring beam. Calibrator J2052+3635 is from the VLBA Calibrator Survey (VCS) bv lKovalev et al.l ([20071) . The other 
calibrators are from the VLA program BZ035. 



a line-free portion of the spectrum. Additionaly, we observed the 
SiO maser emission in two sub-bands centered at Vlsr of - 1 8 
and -5 km s"^ at high spectral resolution (channel spacing of 
97.6563 kHz and bandwidth of 6.25 MHz) with several scans in 
spectral-line mode interspersed among the dual-band continuum 
observations. We also observed H2O maser emission (rest fre- 
quency of 22235.080 MHz) centered at Vlsr of 5 km s"^ with a 
channel spacing of 48.8281 kHz and bandwidth of 6.25 MHz in 
a snapshot mode. 

A typical observing unit for our continuum observations con- 
sisted of a 50 minute NML Cyg scan, followed by a 5 minute 
scan of, alternately, the quasar J2015+3710 and J2012+4628. 
For a spectral-line observation, we observed a calibrator for 
about 3 minutes and NML Cyg for 5 minutes. Absolute flux 
density calibration was established by an observation of 3C286 
(Jl 33 1+3030), assuming a flux density of 1.554 Jy at 43 GHz 
and 2.540 Jy at 22 GHz. Data calibration procedures are de- 
scribed in detail in .Reid & Menten (1997, .2007.) . 

3. Results 

3.1. Maser spectrum and spatial distribution 

Figure [2] shows scalar averaged cross-power spectra of H2O 
maser emission observed using only the inner 5 antennae of the 
VLBA at all five epochs. The H2O maser emission spans a Vlsr 
range from about -25 to 17 km s~\ v v^hich is consistent with 
previ ous observations ( iRichards et all 119961: iNagavama et al ., 
I2008h . One can see that the flux densities of some maser features 
varied considerably from epoch to epoch within about LO yr, 
while the peaks of blue-shifted (Vlsr = -22.17 km s"^) and 
red-shifted (Vlsr = 5.22 km s"^) features appear more stable. 
We confirmed that the brightest and most blue-shifted compo- 
nent appears to exhibit a velocity drift, the LSR velocities of this 
component in 1981, 1993, 2006 and 2008-2009 (this paper) were 
-18, -20, -21.8 and -22.2 km s"^ the resulting acceleration is 
consistent with that derived bv iShintani et al.l (l2008l) . 

Figure [3] shows the spatial distribution of H2O maser emis- 
sion toward NML Cyg relative to the reference maser spot at 
Vlsr= 5.22 km s"^ from observations at five epochs. We con- 
sidered maser spots at diflTerent epochs as being from the same 
feature if their positions in the same spectral channel were co- 
incident within X mas (where x = 3 mas yr"^ x A^ yr, and A^ is 



the time gap between two epochs), corresponding to a linear mo- 
tion of less than 23 km s"^ Selecting the bright maser spots with 
flux density > 0.5 Jy beam"^ and SNR > 10 in each channel, we 
found 21 features that were detected at one or more epochs. 

The total extent of the H2O maser spot distribution is about 
120 mas, whi ch is consisten t with previous studies f rom VLBI 
observations arve i 119961: iNasavama eLaD, l2008l) . Our first 
ep och observation t o ok pla ce about 1 yr later than the last epoch 
in 'Nagayam a et al.l (|2008'), and we found that the most red- 
shifted maser cluster (including features B and E) and the most 
blue shifted maser cluster (including features D and C) in Fig. [3] 
are consistent with features at Vlsr about 6 and -22 km s"^ in 
Nagayama et al. ( 2008), respectively. These two dominant maser 
clusters are about 70 mas apart. This suggests that even though 
the spatial distribution changes in details, the maser features at 
the Vlsr near the distribution of the masers in the two peaks re- 
mains relatively stable. 

Richards et al. (1996) reported an irregular ring of emission 
^ 200 mas across and a pair of outlying features to the north- 
west (NW) and south-ea st (SE), near ly 600 mas a part, based on 
MERLIN observations (lYates & Cohe n. 1994; R ichards et"ail 
[l996) and VLA observations (iJohns'toi et al.L l985). The promi- 
nent maser features in the irregular ring are consistent with the 
features at Vlsr about 5 and -22 km s"^ reported in this paper. 
However, the NW and SE maser features were not detected in 
our VLBA nor in prev ious VLBI observations arve 
Nagayam a et al.L 12008). For our VLBI observations, typical di- 
ameters of the H2O maser spots are ^ 1 mas and the detected lim- 
itation is 10 mJy. This corresponds to a brightness temperature 
lower limit of 4 x 10^ K. However, the brightnes s temper- 
atures for the two outlying features reported by Richa rds et al.l 
(1996) are ^ 10^ K. This might be the reason why the outlying 
features were not detected with VLBI observations. 

Unlike the H2O maser, the SiO maser emission toward NML 
Cyg during our observation was very weak, and phase referenc- 
ing with the maser succeeded only at the last epoch. Because 
we could register H2O and SiO masers to the same background 
sources, the relative position between the reference maser spots 
for two maser species should be accurate to ±1 mas (see Fig.|4l). 
The brightest SiO maser features detected in our observations 
are at Vlsr about -5.70 and 11.03 km s"^ This diflTers from the 
VLBA observatio n of the same v = 1 , / = 1 ^0 SiO maser 
emission by .Boboltz & Marvell (|200Q|) . The extent of the maser 
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Fig. 2. Interferometer (scalar averaged cross-power amplitude over the full duration of the observation) spectra of the H2O masers 
toward NML Cyg obtained with the inner 5 antennae of the VLB A at five epochs. The vertical dash-doted line indicates the maser 
feature at Vlsr of 5.22 km s"^ which served as the phase-reference. The lower panel shows a 100-times blow-up of the H2O spectra 
shown in upper panel. 

(A color version of this figure is available in the online journal.) 



distribution is about 30 mas and is consistent with the diame- 
ter of the ring-like distribution from about 10 yr earlier. ^From 
a VLA observation of v = 0, / = 1 ^0 SiO maser emission, 
iBoboltz & ClaussenI (|2004|) reported a radius of 50 mas for the 
maser distribution, which is consistent with the i nner boundary 
for the dust shell at 2.13 /im (iBlocker et al.Ll200lh . 

3.2. Parallax and proper motion 

We fitted elliptical Gaussian brightness distributions to the im- 
ages of strong maser spots and the extragalactic radio sources 
for all five epochs. The change in position of each maser spot 
relative to each background radio source was modeled by the 
parallax sinusoid in both coordinates (determined by a single 
parameter, the star's parallax) and a linear proper motion in each 
coordinate. The extragalactic source with the best known posi- 
tion, J2052-h3635, was relatively far (3.7°) from NML Cyg, and 
it was only used to determine the absolute position of the maser 
reference spot. We used the other three calibrators, which are 
separated by less than 1.5° from NML Cyg, to determine the 

parallax. 

As mentioned in 'Zhang et aP (1201 2 ^ the apparent motions 
of the maser spots can be complicated by a combination of spec- 
tral blending and changes in intensity. Thus, for parallax fitting, 
one needs to find stable, unblended spots and/or use many maser 



spots to average out these eff'ects. Thus, we first fitted a parallax 
and proper motion to the position off'sets for each H2O maser 
spot relative to each background source separately. In Fig.[5l we 
plot the position of the maser spot at Vlsr = 6.48 km s"^ (the 
first maser spot listed in Table[2l) relative to the three background 
radio sources as an example, with superposed curves represent- 
ing a model that uses a combined estimate for the parallax and 
proper motion as described below. This is one of the most red- 
shifted compact and unblended spots that was detectable at all 
five epochs. 

Table [2l shows the independently estimated parallaxes and 
proper motions for each maser spot and extragalactic source pair. 
While the parallaxes should be identical within measurement un- 
certainties, the proper motions are expected to vary among maser 
spots owing to internal motions of 15 km s"^ (1.9 mas yr"^ at 
a distance of 1.6 kpc). While most of the maser spot parallaxes 
show good internal consistency, the dispersion of parallax esti- 
mates over the entire ensemble is considerably larger than the 
formal errors would suggest. This is caused by residual system- 
atic errors aff'ecting the fits, which originate in the complexity 
and evolution of blended spectral and spatial structure for some 
of the masers. 

Among these independent solutions, we found the parallaxes 
associated with maser spot at Vlsr = -19.64 km s"^ had the 
largest relative uncertainties of ^ 200 yuas. The spots in this fea- 
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Fig. 3. Spatial distribution of the H2O maser features toward NML Cyg from VLB A observations at five epochs. Each maser spot 
is represented by a letter label and di filed circle whose area is proportional to the logarithm of the flux density, using the position 
and flux density from the epoch at which it was first detected. The color bar denotes the Vlsr range from -25 to +25 km s"^. The 
reference maser spot is in feature B, located at (0,0) mas. 
(A color version of this figure is available in the online journal.) 



ture are likely aff'ected by line blended, and we discarded the 
data for this feature. We used all remaining position data in a 
combined solution, which used a single parallax parameter for 
all maser spots and extragalactic source pairs (but allowed for 
diff'erent proper motions for diff'erent maser spots). 

The combined parallax estimate is 0.620 ± 0.047 mas, corre- 
sponding to a distance of 1.61^q|^ kpc, placing NML Cyg pos- 
sibly in a Local (Orion) arm. The quoted uncertainty is the for- 
mal error multiplied by ^[n (where n = 18 is the number of 
maser spots used in the final parallax fit) to allow for the pos- 
sibility of correlated position variations for all the maser spots. 
This could result from small variations in the background source 
or from un-modeled atmospheric delays, b oth of which woul d 
aff'ect the maser spots nearly identically (iReid et al.L l2009al) . 
The un- weighted average absolute proper motion of all selected 
maser spots is /^x = -1.55 ±0.42 mas yr"^ andyu^ = -4.59 + 0.41 
mas yr"\ where j^x = J^a cos S and fiy = /^s- 

Combining the parallax and proper motion measurements 
with the systemic Vlsr of -LO + 2 km s"^ (Kemper et al., 
120031: lEtoka & Diamond! [2004 ) enable us to determine the three- 
dimensional peculiar motion (relative to circular motion around 
the Galactic center) of NML Cyg. Adopting a flat rotation curve 
for the Milky Way with rotation speed of LSR ©o = 239 + 7 
km s~\ distance to the Galactic center of Rq = 8.3 + 0.23 kpc 
(iBrunthaler et al.Ll201 ll) . and the solar motion of ([/0=11.1+^^^, 
V(:)=12.24^^ -|!^, Wp^^ ^'p) km s"^ from Hipparcos measurement 
revised by S chonrich e t al. (2010), we estimate a peculiar motion 
for NML Cyg of (t/,=-2.5 + 3.8, y,=-6.1 + 2.1, W,=-4.8 + 3.3) 
km s"\ where Us, Vg, Wg are directed toward the Galactic cen- 



ter, in the direction of Galactic rotation and toward the North 
Galactic Pole (NGP), respectively. Thus, NML Cyg is orbiting 
the Galaxy slower than expected from a circular orbit for a flat 
rotation curve, similarly t o the aver age of many massive star 
forming regions found by Reid et al. (2009b). The peculiar ve- 
locity components toward the Galactic center and the NGP are 
relatively small with somewhat larger uncertainties. 

3.3. Radio photosphere and circumstellar maser emission 
from VLA observation 

We detected the radio photosphere of NML Cyg at 43 GHz and 
made a uniformly weighted image (see Fig.|7l). We also fitted a 
model of an uniformly bright disk to the wv-data directly using 
the AIPS task OMFIT, obtained a diameter of 44 + 16 mas and 
integrated flux density of 0.82 + 0.11 mJy. We binned the data 
and plotted the real part of the visibility versus baseline length 
in Fig.[6l 

Taking the crude characteristics of the radio photosphere to 
be 7?rp 2 X 7?op, and Trp ~ 0.71 x Tpp, which are inferred from 
the imaging of variable stars by iReid & MentenI {1997) (where 
Rrp and 7?op are the radii of the radio and optical photospheres, 
respectively, and Trp and Top are the efl'ective temperatures of 
the radio and optical photospheres), the semi-diameter of radio 
photosphere, 2 2 mas, is reasonably consistent with that of 16.2 
mas derived bv [Blocker et al.l (l200lh . 

In order to register the SiO masers relative to the radio pho- 
tosphere, we first measured the position of the radio photosphere 
from the broad-band data relative the maser emission from the 
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Fig. 4. Spatial distribution of the H2O (circle ) maser and SiO (square) maser spots toward NML Cyg from the fifth epoch observa- 
tion. Each maser spot is represented by a filled circle or square whose area is proportional to the logarithm of the flux density. The 
position of the central star is indicated by a star (see the text in § 13. 3b . The color bar denotes the Vlsr range from -25 to 25 km s"^ 
of the maser features. The reference H2O maser spot is located at (0,0) mas. 
(A color version of this figure is available in the online journal.) 
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Fig. 6. Fringe visibility vs. baseline length measured for NML 
Cyg at 43 GHz on 2008 December 20 under VLA program 
AZ178. The solid line indicates a uniformly bright circular disk 
model with diameter of 44 mas. 



narrow-band data. Next, we aligned the maser emission from the 
spectral-line data by producing a map of pseudo-continuum data 
by integrating the line data over the entire narrow-band. Because 
the narrow-band and the pseudo-continuum data cover the maser 
emission from the same velocity range, we can align them by 
comparing the positions of peaks in each map. This allows the 
positions of the emission in individual channels to be registered 
to the map of narrow-band data, and then to the radio photo- 
sphere. We estimated a position off'set of the weak continuum 
source (radio photosphere) relative to the reference maser spot 
to be (Ax =-8, A3; = -8) mas with an uncertainty of 3 mas. 



Fig. [7] shows the relative positions of the SiO maser and 43 GHz 
continuum emission used to register the images. 

The relative off'sets of the peaks of SiO maser emission and 
the c entral star agree well w ith observations of Mira variable 
stars (iReid & MentenL |2Q07|) , which show the SiO maser emis- 
sion surrounds the central star in a partial ring with radii of 2 - 
3 /?op. For NML Cyg the Si O radius is also consi stent with the 
dust inner-radius of 45 mas ('Ridgway et al.L [T986b . 

In order to estimate the absolute position of the radio photo- 
sphere (which should be that of the central star), we have to align 
the SiO maser spots measured with the V LA with tho s e from 
the VLB A observations. As discussed in IZhang et aP (l2012h . 
if a single spectral channel contains emission from across the 
source, the VLA position would be falsely interpreted as being 
near the center of the distribution. Furthermore, the VLBA ob- 
serving date of SiO emission is about five months later than that 
of VLA. So, there is no guarantee that the VLA and VLBA maps 
should be identical. However, we found that the two maser clus- 
ters at Vlsr about 11 and -6 km s"^ in the VLBA map are in 
good agreement with those at similar Vlsr values in the VLA 
map, after shifted the VLBA map (Ax = -13, Aj = 26) mas, 
with an uncertainty of 5 mas. Thus, we can register the radio 
photosphere to the VLBA SiO maser spots with a conservative 
uncertainty of 10 mas (see Fig. (8]). 

The absolute position of the reference SiO maser spot can be 
estimated relative to the extragalactic sources, yielding the abso- 
lute stellar position of NML Cyg at epoch 2008.868 of aj2ooo = 
20M6^25?5382 ± 0?0010, ^j2ooo = 40°06'59:'379 ± 0:'015. The 
absolute position is consistent with, but considerably more ac- 
curate than, that from the Two Micron All Sky Survey (2MASS) 
by Skruts kie et al.t ( 2006) (after converting to J2000 coordinates 
and accounting for our measured absolute proper motion); dif- 
ferences of 39 mas in right ascension and 10 mas in declination 
are well within the 170 mas uncertainty of the 2MASS position. 

Fig. [9] shows the H2O maser distribution from our VLA ob- 
servation. Similar to the map as shown in Fig. O the maser dis- 
tribution is dominated by two maser clusters with a separation of 
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Fig. 5. Parallax and proper motion data (markers) and a best-fitting model (dotted lines) that uses a combined estimate for the 
parallax and proper motion for the maser spot at the Vlsr of 6.48 km s"^ Plotted are positions of the maser spot relative to the 
extragalactic radio sources J2044+4005 (circles), J2046+4106 (squares) and J2049+4118 (triangles). Left panel: Eastward (upper 
panel) and northward (lower panel) off'sets versus time. Right panel: Same as the left panel, except the best-fitting proper motion 
has been removed, displaying only the parallax signature. 
(A color version of this figure is available in the online journal.) 



^ 70 mas. However, the total extent of the maser distribution is 
much larger than that of VLB A map. The outlying maser spot to 
the south-east in the VLA map is likely associated with a feature 
reported by .Richards et al. (1996.) . 

4. Discussion 

4.1. Relation to the Cyg region and OB association 

iHab ing et al. (1982) observed a H ii region located northwest of 
NML Cyg, and Morris & Jura (1983) suggested that it repre- 
sented a portion of NML Cyg's circumstellar envelope photoion- 
ized by the luminous and hot stars in the Cyg 0B2 association. 
Cyg OB 2 lies near the center of the X-ray emitting Cygnus X 
super bubble, which is composed of numerous individual H ii 
regions, a number of Wolf-Rayet and 03 stars and several OB 
associations. 

Recently jRygl et all (l2012h reported distances from trigono- 
metric parallax measurements of four 6.7 GHz methanol maser 
sources in the Cygnus X region that are consistent with 1 .40 ± 
0.08 kpc for this region (see Fig. [TT]). Th is is similar to the dis- 
tance of 1 .45 kpc to Cyg 0B2 obtained bv lHansonI (l2003h . based 
on fitting 35 dwarf cluster members with spectral types between 
07.5 and Bl to a 2 x 10^ yr isochrone, and also consistent with 
latest estimated photome tric distance of 1.5 kpc determined by 
iKharchenko et al.l (l2005h . a lthough slightly smaller than previ- 
ous work quoting 1.74 kpc (iMassey & Thomps"on[|l99ll) . 

Given that the Cygnus X region contains huge amounts of 
dust, one might expect that this dust would extinguish the opti- 
cal light from Cyg OB 2. However, little extinction toward Cyg 
OB 2 was found from the NIR extinction map of the region 



(iMotte et al.L l2007l) . We also note that the interstellar medium 
(ISM) absorption signatures from NML Cyg, seen in Herschel 
spectra of the H2O or tho (li.p --> Ip.i) and para ground- state 
(li,o Onn) line (Tev^ier et aLlfin prep .) and low-/ CO tran- 
sitions (iKemper et al .1120031) , suggest that they arise from a cold 
gas component. One possible explanation is that the diff'use gas 
in Cygnus X is located in front of NML Cyg. Taking our dis- 
tance of 1.61^q|^ kpc to NML Cyg, we off'er new observational 
evidence that NML Cyg might be associated with Cyg OB 2. 

As listed in Table [3l our measured absolute proper motion 
of NML Cyg is in good agreement with that of Cyg OB 2, com- 
puted by averaging proper rnotions of the most probable cluster 
members dKharchenko et al.Ll2005l) . NML Cyg's Vlsr of -1 ± 2 
km s"^ is consist with the systemic Vlsr of 7 km s"^ (helio- 
centric velocity -10.3 km s~^), within a dispersion of 6 km s"\ 
for Cyg 0B2 as estimated by Kiminki et atl (l2QQ7l l2008h . For 
comparison. Table [3] also lists parameters for two other red hy- 
pergiants and their associated clusters. We find that the radio 
proper motions of stars and the optical proper motion of their 
associated clusters are in good agreement, providing new obser- 
vational evidence for their association. 

Figure [TT] shows the positions of maser sources in the 
Cygnus X region with accurate parallax measurements and their 
peculiar motions projected on the sky. We find that the methanol 
sources are moving upward (toward increasing Galactic lati- 
tude) while NML Cyg and Cyg OB 2 are moving downward 
from the Galactic plane. This suggest that the space motions in 
the Cygnus X region are much more complicated than expected 
from an expanding Stromgren s phere centered on Cyg OB 2 with 
an angular diameter of 4° dKnodlsedeii l2QQ3h . 
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Fig. 7. SiO maser emissions integrated over Vlsr from -28.86 to -7.14 km s"^ {dashed contour) which is within the same Vlsr 
range as that of the narrow-band data, and Vlsr from -7.14 to 22.93 km s"^ {dotted contour) toward NML Cyg superposed on the 
43 GHz continuum emission {solid contour) with peak position indicated by the star. The maser spots within the same Vlsr range as 
the narrow-band data are marked with circles, with size proportional to the logarithm of their flux densities. The Vlsr of the maser 
spots is color-coded as indicated by the color bar on the right side, which use the same range of Vlsr as other maser spots figures in 
this paper for comparison. All off'sets are relative to the strongest maser emission at Vlsr of -4.66 km s"^ . The relative positions of 
the SiO maser and radio continuum emission is estimated to be accurate to ^ 3 mas. Contour levels for integrated maser emission 
are 0.18 Jy beam"^ x 2", ^ = • • • 5 and for the radio continuum are integer multiples of 0.10 mJy beam"^. The restoring beam for 
the maser {dotted circle) and continuum {solid circle) emissions are indicated in the lower left corner. 
(A color version of this figure is available in the online journal.) 

Table 3. Comparison of astrometric parameters of red hypergiants and its associated clusters 



Star 


Cluster 


Separation 


Distance 






Vlsr 






(°) 


(pc) 


(mas yr~^) 


(mas yr~^) 


(km s-i) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NML Cyg 






1610 + 120 


-1.55 + 0.42 


-4.59 + 0.41 


-01.0 




Cyg 0B2 


2.7 


1500 


-1.22 + 0.20 


-4.72 + 0.24 


+07.0 + 6.0 


VY CMa 






1200 + 120 


-2.80 + 0.20 


+2.60 + 0.20 


+22.0 




NGC2362 


1.3 


1389 


-2.03 + 0.31 


+2.47 + 0.34 


+06.2 + 8.9 


SPer 






2400 + 100 


-0.49 + 0.23 


-1.19 + 0.20 


-38.5 




hPer 


1.5 


2079 


-0.27 + 0.26 


-0.90 + 0.39 


-35.7 + 2.0 




;tPer 


1.5 


2345 


-0.46 + 0.24 


-0.23 + 0.50 


-35.8 + 7.1 



Notes. The first and second columns list the names of a red hypergiant and its associated cluster. The third column lists the separation between 
a star and the center of the cluster. The fourth to seventh columns list distance, proper motion eastward and northward and LSR velocity of red 
hype rgiant or cluster. Par ameters for NML Cyg are from this paper; parame ters for VY CMa are from IZhang et aTl d2012l), except that Vlsr is 
from lMenten et al.l (l2006h : P arameters for S Per are fromlAsaki et al.l (|20TqI) . Parameters for clusters are from rKharchenko et al.l ?2005h . except 
that Vlsr of Cyg 0B2 is from lKiminki et all (l2007lJ2008h . 
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Fig. 8. Comparison of SiO maser distribution from the VLBA observation (squares) on 2009 May 25 and VLA observation (circles) 
on 2008 December 20, after cross-registration (see the text). The sizes of the markers are proportional to the logarithm of their flux 
densities. The Vlsr of the maser spots is color-coded as indicated by the color bar on the right side. The position of the central star 
is indicated a s the s tar. The solid circle and dotted circle denote the stellar diameter of 22 mas from this paper and 16 mas from 
iBlocker et all (l200lh . respectively. 

(A color version of this figure is available in the online journal.) 




East offset (mas) 



Fig. 9. H2O maser emissions distribution toward NML Cyg from VLA observation on 2008 December 20. All off'sets are relative 
to the strongest maser emission at Vlsr of -22.25 km s"^ All the maser spots are marked with circles, with sizes proportional to 
the logarithm of their flux densities. The position of the central star is indicated by a star. 
(A color version of this figure is available in the online journal.) 
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Fig. 10. Evolutionary tracks of 25 and 40 M© non-rotating {dotted lines) and rotating {solid lines) with a velocity of 300 km s"^. 
The model assumes a star with near solar metallicity (Z = 0.02) evolving from the zero-age main sequence to the red supergiant 
and t hen Wolf -Ray et star p hase. The dot-dashed line denotes the empirical upper luminosity boundary (for Teff l ess than 6000 K) 
from lSchuster et al.l (l2006l) . The filled circle labeled "VY CM a" represents the results for that star from.Choi et al.1 (|2Q()8). For NML 
Cy^ t he filled circle labeled "A" corresponds to T^^= 2500 K (iRowan-Robinson & HarrisLll983l:lMonnier et al.Lll997l:IZubko et al.L 
120041) and 2.7 x 10^ L© revised with the distance in this paper; filled circles "B" and "C" are for reff= 2500 K and 3.0 x 10^ L© 
(ISchuster et all. [2009.) and T^^= 3250 K dRMgwav etaL .1986) and 2.7 x 10^ L©. 



4.2. Stellar parameters 

ISchuster et al.l (l2009h determined NML Cyg's minimum bolo- 
metric luminosity, Lboi = (1.04 ± 0.05) x 10^- (J/kpc)^ L©, by 
integrating its SED from 0.5 to 100 jim without correction for 
extinction. This luminosity is in good agreeme nt with the esti - 
mate of Lboi = 1.13 x 10^- ((i/kpc)^ L© by BlockeretaLl (l200lh . 
Adopting our measured distance of 1.61^q|^ kpc, we find NML 
Cyg's luminosity to be (2.7 ± 0.5) x 10^ L©. We note that a 
stellar t emperat ure of 25 CI O K, fi rst derived from a model by 
lRowan- Robinsoi r& HarrisI (1 19831) , gives a good fit to the in- 
frared data (Monnier et al., 1997), as compared to a temperature 
of 3250 K obtained from spectroscopy by Ridgway et al. (1986). 

To examine the evolutionary status of NML Cyg, we com- 
pared i ts position on the H-R d iagram with stellar evolutionary 
tracks (iMeynet & Maedeil 12003 ). using diff'erent Teff and lumi- 
nosity values (see Fig.[TOl). Although NML Cyg's position with 
Teff of 3250 K is more consistent with the evolutionary track 
than those with Teff of 2500 K (compare dots 'C with 'A'), 
we find that NML Cyg's position is closest to the evolutionary 
track with an initial m ass of 25 M©, instead of 40 M© quoted in 
ISchuster et al.l (l2009l) . no matter which luminosity and Teff we 
adopted (see dots 'A, 'B' and 'C in Fig.[TOl). This mass is simi- 
lar to that o f VY CMa, another red hypergaint star, estimated by 
IChoi et afl (|2008). Re-scaling NML Cyg's luminosity (compare 
dots 'A' with 'B') with our distance drops it in the H-R diagram 



well below the empirical upper luminosity boundary (see the red 
dot-dashed line in Fig. [TOj), which is not sensitive to tempera- 
ture for stars < 6000 K. However, NML Cyg's position on the 
H-R diagram, compared with the evolutionary track, constrains 
its age (by adopting the evolutionary time of the closest point to 
NML Cyg on the track) to be 8 Myr, which i s considerably 
greater than that of 2 - 3 Myr for Cyg OB 2 (> Massev et"a]l 
[2001; Hanson, 2003|). Since the observational evidence now con- 
firms that NML Cyg is associated with Cyg 0B2, they are prob- 
ably about the same age. If so, the stellar evolutionary model 
cannot be appropriate for NML Cyg. 

NML Cyg's stellar size of 16.2 mas from Blocker et al.l 
(2001) was derived using the Stepan-Boltzmann law, adopting 
Teff = 2500 K and a distance of 1.74 kpc. Rescaling this stellar 
diameter with our distance of 1.61 kpc gives 15.0 mas. For the 
radio photosphere, we fitted a round-disk model with a diameter 
of 44 ± 16 mas and flux density of 0.8 ± 0.1 mJy using the AIPS 
task OMFIT, and we calculate a (Rayleigh- Jeans) Teff between 
100 K and 6000 K (i.e., ±2cr). Clearly more observations are 
needed to improve this temperature constraint. 

4.3. Maser distribution relative to radio photosphere 

iRichards et al.l (Il996l) proposed a bipolar H2O maser outflow 
and suggested the central star is located close to the prominent 
H2O maser cluster at Vlsr about -22 km s"^ in the irregular 
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Fig. 11. Midcourse Space experiment (MSX) 8 yum image of the Cygnus X region overlaid with the peculiar motions for maser 
sources with parallax and proper motion measurements. T he dotted line indi cates the Galactic plane. The dots mark the water mase r 
(NML Cyg, this paper) and the methanol maser sources ( iRygl et all 1201 2h and the center of Cyg OB 2 (iKharchenko et al.L l2005h . 
with their distances labeled in kpc. The peculiar motions of sources are indi cated with arrows. The approximate size of the Cyg 
0B2 association is denoted by a circle with diameter of 2° (iKnodlsedeil |2QQQ|) . 
(A color version of this figure is available in the online journal.) 



ring, assuming the NW and SE feature lie approximately sym- 
metrically about the st ar. Based on high-resolution HST obser- 
vations, Schuste r et al. |(|2006) find that the circumstellar material 
of NML Cyg shows an "arc-like" shape with a line of symmetry 
and an axis facing towar d Cyg OB 2, which is coincident with 



the H2O maser (Richards et al.Ul996| ) and ground state SiO dis- 
tribution (Boboltz & Claussen, 2004) around NML Cyg. They 
assumed the star is near the peak intensity in the optical images 
and coincident with the position of the strongest H2O maser fea- 
ture, then the maser shows an one-sided distribution within NML 
Cyg's nebula. This led to the suggestions that the masers are pro- 
tected by the star's envelope from Cyg 0B2's radiation and the 
arc-like shape of circumstellar material is the result of the inter- 
action between the molecular outflow from NML Cyg and th e 
near-UV continu um flux from Cyg 0B2 ( Schust er et alll2009l) . 
iNagayama et al. l ( 2008) presented three-dimensional kinematic 
data for H2O maser features from multi-epoch VLBI observa- 
tions and suggested an expanding outflow with the central star at 
the dynamical center. 



However, after registering the radio photosphere to the SiO 
maser and then to the H2O maser, we find that both the SiO 
maser and H2O maser are highly asymmetrically distributed 
around the NML Cyg's central star, which is located at south- 
west side of both the H2O and SiO maser distribution, and there 
are H2O and SiO masers located north-west of th e central star. 
This is inconsistent with the model suggested by Schu ster et al.l 
(2006, 2009) based on the one-side H2O maser distribution, 
since there are no masers located north-west of the central star. 
And the central star is not near the center of the two promi- 
nent maser clus t ers at Vlsr of about -22 and 5 km s"^ as 
Nagavam a et in (l2QQ8h suggested, but is closer to the maser 
cluster at Vlsr of about -22 km s~\ w hich is coincident with 
that suggested bv lRichardsetaD(ll996h . Considering that kine- 
matic models for maser features are highly dependent on maser 
positions and motions respective to the central star, our direct 
determination of positions relative to the central star presents a 
crucial constraint for understanding the properties of the circum- 
stellar material surrounding NML Cyg. 
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5. Conclusions 

We have measured the trigonometric parallax and proper mo- 
tion of NML Cyg from multiple epoch VLB A observations of 
the circumstellar H2O and SiO masers. Both the distance and 
proper motion of NML Cyg are consistent with those of the Cyg 
0B2 association within their joint uncertainties, confirming that 
NML Cyg is associated with Cyg 0B2. We revised the stellar 
luminosity using the accurate distance from a trigonometric par- 
allax measurement, suggesting that the location of NML Cyg on 
the HR diagram is consistent with the evolutionary track of an 
evolved star with an initial mass of 25 M©. However, the derived 
age of 8 Myr for NML Cyg from the evolutionary track is much 
larger than that 2-3 Myr of Cyg OB 2, suggesting that the evolu- 
tionary track might not be appropriate to NML Cyg. 

After registering the radio photosphere to the H2O and SiO 
maser, we find the central star is located close to a prominent 
H2O maser cluster towards the north-west edge of the maser dis- 
tribution. This is inconsistent with the model that the masers are 
off to one- side of the star and also the model for an expanding 
outflow with the star in the middle of two prominent H2O maser 
clusters as previously suggested. 
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Table 2. Parallax and proper motion fits 



B. Zhang et al.: Distance and size of NML Cyg 



Source 



Region 



^LSR 

(kms-1) 



Parallax 

(mas) 



(mas yr"^) 



(mas yr ^) 



Ax 
(mas) 



Ay 
(mas) 



J2044+4005 


B 


6.48 


0.68 


+ 


0.10 


-1.66 


+ 


0.58 


-4.53 


+ 


0.19 


-0.44 


+ 


0.19 


2.19 


+ 


0.06 




B 


6.06 


0.72 


± 


0.10 


-1.70 


± 


0.52 


-4.33 


± 


0.19 


-0.37 


± 


0.17 


2.09 


± 


0.05 




B 


5.64 


0.70 


+ 


0.09 


-1.59 


+ 


0.48 


-4.47 


+ 


0.18 


-0.50 


+ 


0.16 


2.18 


+ 


0.05 




B 


5.22 


0.70 


+ 


0.08 


-1.68 


+ 


0.44 


-4.55 


+ 


0.15 


-0.63 


+ 


0.15 


2.29 


+ 


0.04 




B 


4.80 


0.69 


± 


0.06 


-1.31 


± 


0.47 


-4.73 


± 


0.12 


-0.90 


± 


0.16 


2.55 


± 


0.03 




B 


4.37 


0.55 


+ 


0.06 


-1.85 


+ 


0.18 


-5.50 


+ 


0.15 


-0.66 


+ 


0.06 


3.10 


+ 


0.04 




B 


3.95 


0.73 


+ 


0.07 


-1.08 


+ 


0.52 


-4.59 


+ 


0.14 


-0.25 


+ 


0.17 


3.80 


+ 


0.04 




F 


-18.80 


0.69 


± 


0.10 


-1.81 


± 


0.43 


-5.25 


± 


0.20 


-85.02 


± 


0.14 


0.71 


± 


0.06 




F 


-19.22 


0.65 


+ 


0.12 


-1.90 


+ 


0.39 


-5.43 


+ 


0.25 


-84.93 


+ 


0.13 


0.94 


+ 


0.08 




F 


-19.64 


0.88 


+ 


0.21 


-3.34 


+ 


0.93 


-5.01 


+ 


0.41 


-83.88 


+ 


0.31 


6.87 


+ 


0.12 




G 


-20.06 


0.59 


± 


0.12 


-1.28 


± 


0.33 


-4.37 


± 


0.32 


-83.05 


± 


0.11 


6.98 


± 


0.10 




A 


-20.48 


0.66 


+ 


0.09 


-1.99 


+ 


0.36 


-4.85 


+ 


0.18 


-83.15 


+ 


0.12 


6.99 


+ 


0.05 




A 


-20.90 


0.73 


± 


0.10 


-2.25 


± 


0.45 


-4.98 


± 


0.19 


-83.20 


± 


0.15 


6.96 


± 


0.05 




A 


-21.32 


0.66 


± 


0.09 


-1.24 


± 


0.33 


-5.06 


± 


0.20 


-76.67 


± 


0.11 


5.59 


± 


0.06 




A 


-21.75 


0.66 


+ 


0.10 


-1.31 


+ 


0.32 


-4.99 


+ 


0.22 


-76.78 


+ 


0.11 


5.63 


+ 


0.07 




A 


-22.17 


0.75 


± 


0.07 


-2.23 


± 


0.62 


-5.15 


± 


0.14 


-84.58 


± 


0.20 


4.52 


± 


0.04 




A 


-22.59 


0.75 


± 


0.08 


-2.06 


± 


0.54 


-5.15 


± 


0.14 


-84.71 


± 


0.18 


4.50 


± 


0.04 




A 


-23.01 


0.70 


+ 


0.09 


-2.12 


+ 


0.42 


-5.17 


+ 


0.17 


-84.77 


+ 


0.14 


4.46 


+ 


0.05 




R 


-24.27 


0.71 


± 


0.06 


-0.63 


± 


0.53 


-3.88 


± 


0.12 


-41.28 


± 


0.18 


33.57 


± 


0.03 


J2046+4106 


B 


6.48 


0.63 


+ 


0.05 


-1.56 


+ 


0.13 


-4.29 


+ 


0.13 


-0.23 


+ 


0.05 


2.73 


+ 


0.04 




B 


6.06 


0.62 


± 


0.01 


-1.59 


± 


0.02 


-4.12 


± 


0.16 


-0.16 


± 


0.01 


2.64 


± 


0.05 




B 


5.64 


0.64 


± 


0.01 


-1.50 


± 


0.02 


-4.24 


± 


0.15 


-0.29 


± 


0.01 


2.72 


± 


0.05 




B 


5.22 


0.69 


+ 


0.01 


-1.60 


+ 


0.02 


-4.27 


+ 


0.17 


-0.42 


+ 


0.01 


2.83 


+ 


0.06 




B 


4.80 


0.61 


± 


0.01 


-1.21 


± 


0.02 


-4.53 


± 


0.22 


-0.68 


± 


0.01 


3.10 


± 


0.07 




B 


4.37 


0.57 


± 


0.09 


-1.77 


± 


0.26 


-5.20 


± 


0.22 


-0.46 


± 


0.09 


3.64 


± 


0.07 




B 


3.95 


0.56 


+ 


0.03 


-0.96 


+ 


0.06 


-4.46 


+ 


0.25 


-0.02 


+ 


0.02 


4.36 


+ 


0.08 




F 


-18.80 


0.62 


± 


0.03 


-1.72 


± 


0.08 


-5.03 


± 


0.18 


-84.81 


± 


0.03 


1.26 


± 


0.06 




F 


-19.22 


0.64 


± 


0.03 


-1.82 


± 


0.07 


-5.14 


± 


0.07 


-84.73 


± 


0.02 


1.49 


± 


0.02 




F 


-19.64 


0.81 


+ 


0.20 


-3.24 


+ 


0.50 


-4.77 


+ 


0.59 


-83.67 


+ 


0.17 


7.42 


+ 


0.19 




G 


-20.06 


0.61 


± 


0.10 


-1.21 


± 


0.28 


-4.04 


± 


0.26 


-82.85 


± 


0.09 


7.52 


± 


0.08 




A 


-20.48 


0.65 


± 


0.04 


-1.91 


± 


0.11 


-4.57 


± 


0.15 


-82.94 


± 


0.04 


7.54 


± 


0.05 




A 


-20.90 


0.67 


+ 


0.02 


-2.15 


+ 


0.04 


-4.75 


+ 


0.17 


-82.99 


+ 


0.02 


7.51 


+ 


0.06 




A 


-21.32 


0.64 


± 


0.04 


-1.16 


± 


0.10 


-4.78 


± 


0.21 


-76.47 


± 


0.03 


6.14 


± 


0.07 




A 


-21.75 


0.64 


± 


0.05 


-1.23 


± 


0.11 


-4.71 


+ 


0.19 


-76.57 


± 


0.04 


6.17 


+ 


0.06 




A 


-22.17 


0.53 


± 


0.01 


-2.10 


± 


0.02 


-5.07 


± 


0.20 


-84.34 


± 


0.01 


5.08 


± 


0.06 




A 


-22.59 


0.58 


± 


0.01 


-1.95 


± 


0.02 


-5.02 


± 


0.18 


-84.48 


± 


0.01 


5.05 


± 


0.06 




A 


-23.01 


0.63 


± 


0.09 


-2.03 


± 


0.38 


-4.94 


± 


0.17 


-84.56 


± 


0.13 


5.01 


± 


0.05 




R 


-24.27 


0.54 


+ 


0.02 


-0.52 


+ 


0.05 


-3.75 


+ 


0.28 


-41.05 


+ 


0.02 


34.13 


+ 


0.09 


J2049+4118 


B 


6.48 


0.62 


± 


0.02 


-1.48 


± 


0.07 


-4.03 


± 


0.05 


0.13 


+ 


0.02 


2.15 


± 


0.02 




B 


6.06 


0.63 


+ 


0.03 


-1.52 


+ 


0.08 


-3.82 


+ 


0.13 


0.20 


+ 


0.03 


2.05 


+ 


0.04 




B 


5.64 


0.64 


± 


0.04 


-1.42 


± 


0.10 


-3.95 


± 


0.12 


0.07 


± 


0.03 


2.13 


± 


0.04 




B 


5.22 


0.68 


± 


0.05 


-1.52 


± 


0.11 


-3.99 


± 


0.16 


-0.06 


± 


0.04 


2.24 


± 


0.05 




B 


4.80 


0.62 


+ 


0.04 


-1.14 


+ 


0.09 


-4.23 


+ 


0.18 


-0.33 


+ 


0.03 


2.51 


+ 


0.06 




B 


4.37 


0.52 


± 


0.07 


-1.68 


± 


0.35 


-4.95 


± 


0.14 


-0.09 


± 


0.11 


3.06 


± 


0.04 




B 


3.95 


0.59 


± 


0.06 


-0.89 


± 


0.15 


-4.14 


± 


0.19 


0.33 


± 


0.05 


3.77 


± 


0.06 




F 


-18.80 


0.64 


+ 


0.06 


-1.64 


+ 


0.17 


-4.72 


+ 


0.13 


-84.45 


+ 


0.06 


0.67 


+ 


0.04 




F 


-19.22 


0.62 


± 


0.01 


-1.73 


± 


0.13 


-4.86 


± 


0.03 


-84.37 


± 


0.04 


0.90 


± 


0.01 




F 


-19.64 


0.81 


+ 


0.17 


-3.17 


+ 


0.42 


-4.47 


+ 


0.53 


-83.31 


+ 


0.14 


6.83 


+ 


0.17 




G 


-20.06 


0.62 


+ 


0.10 


-1.13 


+ 


0.27 


-3.74 


+ 


0.26 


-82.49 


+ 


0.09 


6.93 


+ 


0.08 




A 


-20.48 


0.62 


± 


0.04 


-1.83 


± 


0.16 


-4.32 


± 


0.08 


-82.58 


± 


0.05 


6.95 


± 


0.02 




A 


-20.90 


0.67 


± 


0.04 


-2.08 


+ 


0.10 


-4.45 


+ 


0.12 


-82.63 


+ 


0.03 


6.92 


± 


0.04 




A 


-21.32 


0.64 


+ 


0.06 


-1.08 


+ 


0.15 


-4.49 


+ 


0.15 


-76.11 


± 


0.05 


5.55 


+ 


0.05 




A 


-21.75 


0.65 


± 


0.06 


-1.15 


± 


0.17 


-4.41 


± 


0.14 


-76.21 


± 


0.06 


5.58 


± 


0.04 




A 


-22.17 


0.55 


+ 


0.05 


-2.03 


+ 


0.11 


-4.76 


+ 


0.21 


-83.99 


+ 


0.04 


4.49 


+ 


0.07 




A 


-22.59 


0.60 


+ 


0.05 


-1.87 


+ 


0.10 


-4.70 


+ 


0.18 


-84.13 


+ 


0.04 


4.46 


+ 


0.06 




A 


-23.01 


0.65 


+ 


0.05 


-1.95 


± 


0.47 


-4.64 


± 


0.09 


-84.20 


± 


0.15 


4.42 


± 


0.03 




R 


-24.27 


0.56 


± 


0.06 


-0.44 


± 


0.13 


-3.44 


± 


0.22 


-40.70 


± 


0.05 


33.54 


± 


0.07 
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Table 2. Continuted. 



B. Zhang et al.: Distance and size of NML Cyg 



Source 



Region 



V'lsr 
(kms-1) 



Parallax 
(mas) 



(mas yr"^) 



(mas yr ^) 



Ax 
(mas) 



(mas) 



Combined 



B 


6.48 


-1.56 + 0.13 


-4.29 + 0.11 


B 


6.06 


-1.59 + 0.13 


-4.11 + 0.11 


B 


5.64 


-1.49 + 0.13 


-4.24 + 0.11 


B 


5.22 


-1.58 + 0.13 


-4.33 + 0.11 


B 


4.80 


-1.21 + 0.13 


-4.50 + 0.11 


B 


4.37 


-1.78 + 0.13 


-5.14 + 0.11 


B 


3.95 


-0.97 + 0.13 


-4.39 + 0.11 


F 


-18.80 


-1.71 + 0.13 


-5.02 + 0.11 


F 


-19.22 


-1.81 + 0.13 


-5.14 + 0.11 


G 


-20.06 


-1.21 + 0.13 


-4.02 + 0.11 


A 


-20.48 


-1.90 + 0.13 


-4.59 + 0.11 


A 


-20.90 


-2.14 + 0.13 


-4.78 + 0.11 


A 


-21.32 


-1.15 + 0.13 


-4.79 + 0.11 


A 


-21.75 


-1.22 + 0.13 


-4.71+0.11 


A 


-22.17 


-2.11 + 0.13 


-4.98 + 0.11 


A 


-22.59 


-1.95 + 0.13 


-4.97 + 0.11 


A 


-23.01 


-2.03 + 0.13 


-4.94 + 0.11 


R 


-24.27 


-0.53 + 0.13 


-3.66 + 0.11 



0.620 + 0.047 



Notes. Absolute proper motions are defined as fi^ - f^a cos 6 and fiy -115. 
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